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We investigate the impact of rotational diffusion on the electrodynamic coupling of fluorescent 
dye molecules (oscillating electric dipoles) to a tunable planar metallic nanocavity Fast rotational 
diffusion of the molecules leads to a rapidly fluctuating mode density of the electromagnetic field 
along the molecules' dipole axis, which significantly changes their coupling to the field as compared 
to the opposite limit of fixed dipole orientation. We derive a theoretical treatment of the problem 
and present experimental results for rhodamine 6G molecules in cavities filled with low and high 
viscosity liquids. The derived theory and presented experimental method is a powerful tool for 
determining absolute quantum yield values of fiuorescence. 
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Introduction. — Fluorescing molecules located close to 
a metal surface (at sub-wavelength distance) or inside a 
metal nano-cavity, dramatically change their fluorescence 
emission properties such as fluorescence lifetime, fluores- 
cence quantum yield, emission spectrum, or angular dis- 
tribution of radiation [IHlj • This is due to the change lo- 
cal density of modes of the electromagnetic field caused 
by the presence of the metal surfaces [S]. Although a 
large amount of studies have dealt with the investigation 
of this effect, they all have considered fixed dipole orien- 
tations of the emitting molecules, so that each molecule 
exhibits a temporally constant mode density during its 
de-excitation from the excited to the ground state. How- 
ever, when molecules arc dissolved in a solvent such as 
water, their rotational diffusion leads to rapid changes 
of dipole orientation even on the time-scale of the av- 
erage excited state lifetime. We will show here that 
this dramatically influences the coupling of the molecules 
to the local, strongly orientation-dependent density of 
modes and the resulting excited state lifetime. This is 
enormously important for applications of tunable nano- 
cavities for fluorescence quantum yield measurements. 

Theory. — Let us consider an ensemble of molecules 
within a planar nano-cavity, which had been excited by 
a short laser pulse into their excited state. Due to the 
electrodynamic coupling to the cavity, these molecules 
will exhibit an emission rate K that depends on their 
vertical position within the cavity, and on the angle 9 
between their emission dipole axis and the vertical. In 
what follows, we assume that the excited state lifetime 
is so short that one can neglect any translational diffu- 
sion of a molecule within the cavity. However, this is 
in general not the case for its rotational diffusion time 
which can be on the same order as the excited state life- 



time. Then, for a given position within the cavity, the 
probability density p{9, t) to find a molecule still in its 
excited state at time t with orientation angle 9 obeys the 
following evolution equation 



dp(9,t) 
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Bpi9,t) 



K{9)p{e,t) (1) 



dt sin 961 89 

where the first term on the right hand side is the ro- 
tational diffusion operator [B] multiplied with rotational 
diffusion coefficient D, and the second term accounts for 
de-excitation. For the sake of simplicity, we omit any 
explicit indication of the position dependence of the in- 
volved variables. The emission rate K itself is given by 
a weighted average of the wavelength dependent rates 
k{9,\), 



K{9) = {k{9,X))x = 



Jk{9,X)Fo{X)dX 



(2) 



jFo{X)dX 

where Fo{X) is the free-space emission spectrum of the 
molecules as a function of wavelength A. For a planar 
cavity, the rates k themselves can be decomposed into 



k{9, A) = k_i_ (A) cos^ 9 + kn (A) sin^ 



(3) 



where the fcj^ || (A) are the rates for a vertically and a 
horizontally oriented emitter, respectively. Within the 
semi-classical theory of dipole emission [71, these rates 
are given by 
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FIG. 1: Computational results for the decadic logarithm of 
the relative lifetime, Ig So/S^{X), as a function of wavelength 
A and vertical position z within a nano-cavity comprising of a 
silver layer of 35 nm (bottom mirror of cavity) and one of 85 
nm (top mirror of cavity). The cavity is filled with water. The 
top panel shows the result for a vertical dipole orientation, the 
bottom panel for a horizontal orientation. The cavity size, i.e. 
distance between silver surfaces, is 173 nm corresponding to a 
white- light majdmum of transmission of the cavity of 601 nm. 
The white line in the bottom panel divides the region where 
the ratio So/Sfj,{X) is larger than one from that where it is 
smaller than one. 



where the index /i is either _L or ||, and where knr and 
krad are the free-space non-radiative and radiative tran- 
sition rates, respectively, tq is the free-space excited state 
lifetime, $ is the intrinsic quantum yield of fluorescence, 
S)j,{X) are the wavelength-dependent emission rates of an 
oscillating electric dipole with orientation /i within the 
cavity, and Sq is the free-space emission rate, which is 
independent on orientation and wavelength (thus neglect- 
ing optical dispersion of the solvent). The emission rates 
Sfj,{X) are calculated in a semi-classical way by firstly us- 
ing a plane wave representation of the electromagnetic 
field of an emitting electric dipole of given orientation 
(and position) [5]; secondly calculating the interaction of 
each plane wave component with the cavity; and finally 
finding the emission rate as the integral of the Poynting 
vector of the total field over two surfaces sandwiching 
the emitter on both sides. An exemplary result for such 
a calculation is shown in Fig. [T] 

The initial distribution p(^, t = 0) right after excitation 
is defined by the polarization and intensity of the focused 
excitation light. These can be found by again expanding 
the electromagnetic field of the focused laser beam into 
a plane wave representation [HI HO] , and calculating the 
interaction of each plane wave with the cavity |11| . If 
one denotes the horizontal and vertical components of 
the excitation intensity at the position of the molecules 
by /|| and Ij_, respectively, then p{9,t = 0) is given by 



p{B,t = {)) = 



3(/_L cos^ e + /,, sin^ 



(5) 



2(/_L + 2/||) 
Computational results for I±^ and In are shown in Fig. [21 
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FIG. 2: Computational results for the normalized excita- 
tion intensity distribution within the same nano-cavity as de- 
scribed in Fig. [l| For better visualization of the cavity's ge- 
ometry, the figure shows also the silver layers as gray-shaded 
areas. The left side of the figure shows the excitation in- 
tensity for horizontally oriented molecules, the right side for 
vertically oriented molecules. 



for the same cavity geometry as in Fig. [T] 

Next, the solution to Eq. ([l]) can be found by expand- 
ing p{9, t) into a series of Legendre polynomials ^^(cos 9): 



p(0,t)=^a,(i)P,(cos0) 



(6) 
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where the ai{t) denote time-dependent expansion coeffi- 
cients. Inserting that into Eq. ([I]) yields an infinite set of 
ordinary differential equations for the ai{t), 



daijt) 
dt 



= -Dii£ + l)ae{t) - ^ Mwa^, (t) (7) 



where the transition matrix Mu' is defined by the inte- 
grals 
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with the abbreviations Kj_^\\ = (^_l,|i(A))a, and AA' — 
K± — Kn . By carrying out the integration, one finds that 
the only non-vanishing components of Mui are given by 
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From the initial condition, Eq. ([s]), one finds that the 
only non-vanishing initial values ai are 



ao{t = 0) 
a2it = 0) 



I± - I\\ 

h +2/,, 



(10) 



Although Eq. ([t]) represents an infinite set of differential 
equations, it occurs that for our experimental conditions 



(see below) a truncation of the series expansion of Eq. ^ 
at a maximum £max = 10 yields an accurate solution to 
the problem that does not change when further increasing 
this truncation value. 

It remains to find an expression for the observable flu- 
orescence emission. This is given by the integral 



Fit) 



dr 



sin0p{0,t){k{e,X)u{0,X))x (11) 



where u{0, A) is the orientation and wavelength depen- 
dent fluorescence detection efficiency, {)\ denotes inte- 
gration over all wavelengths, and the first integration ex- 
tends over the whole inner space of the cavity. Due to 
the rapid fall-off of the excitation intensity when moving 
a few micrometers away from the center of the focused 
laser beam, the integration over space can be cut off ac- 
cordingly. Similarly to the emission rate, the detection 
efficiency can be represented by 



u{9, A) = u±{X) cos^ 9 + Mji (A) sin^ 



(12) 



with u± and uy being the detection efficiencies for a verti- 
cally and horizontally oriented emitter. The most signif- 
icant cause which makes these detection efficiencies dif- 
ferent is the strongly orientation-dependent angular dis- 
tribution of radiation of the emitters which is collected 
differently by the detection optics with finite aperture. 
The detection efficiencies are calculated again via a plane 
wave representation of the emitted electromagnetic field, 
for details see [121 E] . It should be noted that the detec- 
tion efficiency goes down to zero when approaching the 
silver mirrors so that only fluorescence from molecules 
at least a few nanometers away from the cavity surfaces 
contributes to the detected signal. 

When inserting the expansion (IbI) into Eq. (11) and 
integrating over 9 one finds that only the amplitudes a^ 
with £ e {0, 2, 4} contribute to the final result. 



F{t) ^ dr [Coaoit) + C2a2{t) + Cia^it)] (13) 

while the constant factors C^ are given by 



Co — (3fc_LU_L + 2fc||M_L + 2fc_LM|| + 8fc||M||)A 
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C2 =— 3(6fc_LM_L + h\Ul_ + fc_LU|| - 8fc||U||>A (14) 

105 
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C4 =77T-(fcj_iij_ - k\\u±_ - k±u\\ + k\\u\\)x 
315 

Finally, the observable mean fluorescence lifetime (t) is 
found as 



Experiment. — A homemade nano-cavity consists of two 
silver mirrors with sub-wavelength spacing. The bottom 
silver mirror (35 nm thick) was prepared by vapor deposi- 
tion onto commercially available and cleaned microscope 
glass coverslides (thickness 170 ^m) using an electron 
beam source (Laybold Univex 350) under high-vacuum 
conditions (~10~^ mbar). The top silver layer (85 nm 
thick) was prepared by vapor deposition of silver onto the 
surface of a plan-convex lens (focal length of 150 mm) un- 
der the same conditions. Film thickness was monitored 
during vapor deposition using an oscillating quartz unit 
and verified by atomic force microscopy. The complex- 
valued wavelength-dependent dielectric constants of the 
silver films were determined by ellipsometry (nanofilm 
ep3se, Accurion GmbH, Gottingen) and subsequently 
used for all theoretical calculations. The spherical shape 
of the upper mirror allowed us to reversibly tune the cav- 
ity length by retracting from or approaching to the cavity 
center. It should be noted that within the focal spot of 
the microscope objective lens the cavity can be consid- 
ered as a plane-parallel resonator [T3]. For the lifetime 
measurements, a droplet of a micromolar solution of rho- 
damine 6G molecules in water or glycerol was embedded 
between the cavity mirrors. The cavity length was deter- 
mined by measuring the white light transmission spec- 
trum [TH [TS] using a spectrograph (Andor SR 303i) and 
a GGD camera (Andor iXon DU897 BV), and by fitting 
the spectra with a standard Fresnel model of transmis- 
sion through a stack of plan-parallel layers, where the 
cavity length (distance between silver mirrors) was the 
only free fit parameter. 

Fluorescence lifetime measurements were performed 
with a home-built confocal microscope equipped with an 
objective lens of high numerical aperture (UPLSAPO, 
60x, N.A. = 1.2 water immersion, Olympus). A white- 
light laser system (Fianium SG400-4-80) with a tunable 
mter (AOTFnC-400.650-TN) served as the excitation 
source (Aoxc = 488 nm). The light was reflected by a 
dichroic mirror (Semrock BrightLine FF484-FDi01) to- 
wards the objective, and back-scattered excitation light 
was blocked with a long pass filter (Semrock EdgeBasic 
BLP01-488R). Gollected fluorescence was focused onto 
the active area of an avalanche photo diode (PicoQuant 
r-SPAD). Data acquisition was accomplished with a mul- 
tichannel picosecond event timer (PicoQuant HydraHarp 
400). Photon arrival times were histogrammed (bin 
width of 50 ps) for obtaining fluorescence decay curves, 
and all curves were recorded until reaching 10^ counts 
at the maximum . Finally, the fluorescence decay curves 
were fitted with a multi-exponential decay model, from 
which the average excited state lifetime was calculated 



according to Eq. ( 15 1 



(T) = ; dtF{t)t 



dtF{t) 



(15) 



Fig. [3] shows the result of the measured average fluo- 
rescence lifetime of rhodamine 6G in water (blue dots) 
and glycerin (red dots) within the nano-cavity as a func- 
tion of maximum transmission wavelength (which is lin- 
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FIG. 3: Fluorescence lifetime of rhodamine 6G in water and 
in glycerol as a function of the maximum transmission wave- 
length of the cavity (i.e. cavity size). Dots are experimental 
data, solid lines theoretical curves, and dotted lines indicate 
the free-space lifetime values. 



beam will predominantly excite molecules with horizon- 
tal orientation (see Fig. [2]) , for which the emission rate 
can be lower than the free-space rate. If the molecules 
do not rotate, one can thus observe, for specific cavity 
size values, average lifetime values which are longer than 
the free-space lifetime. However, if molecular rotation is 
much faster than the average excited state lifetime, than 
the emission rate will be dominated by that for verti- 
cally oriented molecules (which is much faster than that 
for horizontally oriented ones, see Fig. [I]) and will always 
result in average lifetime values smaller then the free- 
space lifetime. Finally, it should be emphasized that the 
excellent agreement between theoretical model and ex- 
perimental results offer the fascinating possibility to use 
lifetime measurements on dye solutions in tunable nano- 
cavities for simple and direct determination of the fluo- 
rescence quantum yield, a quantity which is notoriously 
difficult to determine by other methods [15]. 

Financial support by the Deutsche Forschungsgemein- 
schaft is gratefully acknowledged (SFB 937, project A5). 



early proportional to cavity length). Both curves show 
a strong decrease of the lifetime values with increasing 
cavity length. The solid lines represent fits of the theo- 
retical model to the experimental values, where the only 
fit parameters have been the free space lifetime tq, the 
fluorescence quantum yield $, and the rotational diffu- 
sion time td = 1/61?. For the water solution, the fit 
values are tq = 4.1 ns, $ = 0.93, and td < 50 ps, 
whereas for the glycerol solution they are tq ~ 3.6 ns, 
$ = 0.99, and t^ > 100 ns (indicating that Z? w on 
the time scale of the fluorescence lifetime). The fluo- 
rescence lifetime and fluorescence quantum yield values 
are in excellent agreement with published values, see [16] 
and citations therein. The large fit value of the rota- 
tional diffusion time for rhodamine in glycerol, which 
is by nearly two orders of magnitude larger than the 
fluorescence lifetime, indicates that rotational diffusion 
is practically frozen during de-excitation of the excited 
molecules, which is similar to the limiting case of fixed 
dipole orientations. Contrary, the fitted rotational diffu- 
sion value in water is significantly shorter than the life- 
time, indicating a situation where the emitters perceive 
an environment with a rapidly fiuctuating mode density 
of the electromagnetic field. Both situations, rapid and 
slow rotational diffusion, lead not only to quantitatively 
different results for the dependence of average lifetime 
on cavity size as seen in Fig. |3l but also to qualitatively 
different behavior: While for slow rotators, the average 
lifetime can exceed, for specific cavity size values, the free 
space lifetime (dotted lines in Fig.lsl, the average lifetime 
for rapidly rotating molecules will always be smaller than 
the free-space lifetime. The reason for that can be under- 
stood when inspecting Figs, [l] and [2j The focused laser 
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